CONVEX AND LINEAR MODELS OF NP-PROBLEMS 



SERGEY GUBIN 



Abstract. Reducing the NP-problems to the convex/hnear analysis on the 
BirkhofT polytope. 



Introduction 

Since the classical works of J. Edmonds ^ and others], linear modeling became 
a common technique in combinatorial optimization [51 [9l [131 [Ml [151 [TBI and others]. 
Often, the linear models are expressed with some constrains on the incidence vector. 
The major benefit of this approach is the symmetry of the resulting model: the 
resulting equations are an invariant under relabeling. The major disadvantage of 
the approach is difficulty to express the constrains explicitly due to their size and 
structural complexity ^ [9l [131 and others] . 

In this work, the Subgraph Isomorphism Problem [3l [H [7l [11] is taken as a 
basic NP-problem. The adjacency and incidence matrices are used to express the 
linear and convex models explicitly. In such asymmetric models, the unknown is a 
relabeling. The relabeling is presented with an unknown permutation matrix. That 
reduces the NP-problems to the linear/convex analysis on the Birkhoff polytope [T]. 

1. Adjacency matrix models 

Let's take the Subgraph Isomorphism Problem [3l|4l[7]: whether a given multi 
digraph g contains a subgraph which is isomorphic with another given multi digraph 
s. That is a NP-complete problem. 

Let n and m be powers of vertex-sets of g and s, appropriately. Based on a 
node labeling/enumeration, let's construct adjacency matrices of these digraphs - 
matrices G and S, appropriately. In terms of these matrices, the problem is a 
compatibility problem for the following quadratic system where the unknown is 
permutation matrix X — {xij)nxn'- 

(1.1) PjnnX^ GX > S 

Here, matrix Pmn is a truncation: 

- where matrix Um is the union matrix of size m x m. 

Permutation matrix X presents all n\ possible ways to label vertices of (multi) 
digraph g. Compatibility of svstem fl . 1 1 means that there is such a vertex labeling 
of digraph g that the appropriate adjacency matrix of that digraph will "cover" 
the adjacency matrix of (multi) digraph s. Such a labeling of g would make the 
positive solution for instance (g, s) the self evident. Obviously, if the system is 
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incompatible, then instance {g,s) has solution "NO". 

Model fTTTl is asymmetric [9] because relabeling of g rotates solutions of system [TTT] 
over all permutation matrices. The matrices are vertices of the Birkhoff polytope. 

Let's call matrices G and S the problem's instance and pattern, appropriately. 
The matrix pairs (or more precisely, the pair's conjugacy classes over the group of 
permutation matrices) parametrize the whole NP zoo. Let us illustrate that with 
several examples [2l|3l|tl[2l[IIl[l4l[15l[l6]: 

(Sub)GI: Instance and pattern are adjacency matrices of given (multi di-) 
graphs. In case of the GI problem: P„j„ = [/„ and sign ">" is replaced 
with sign "=" . 

Clique: Instance is an adjacency matrix of a given (multi di-) graph. Pattern 
is a matrix whose diagonal elements are and remaining elements are 1: 



S = (1). 



1\ 



/ 



(1.2) 



HC: Instance is adjacency matrix of a given (multi di-) graph. Pattern is any 
circular permutation matrix, for example: 

1 \ 





s 



( 

1 





1 0/ 



V 

HP: The same as the above except one 1 is poked out: 

/ ... \ 



\ 



/ 

/ n X n 



Matching: In this case, m is an even number. Pattern S can be, for example, 
the following matrix: 



S 



( 


1 





























1 














V : 









- where 1 and are altering on the over-diagonal. 
Perfect Matching: The same as the above except m = n. 
3-SAT: Let / be a 3-SAT instance: 



f^f\LaV La V La. 
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Let's arbitrarily enumerate strings in the truth tables of the given clauses. 
By definition, two strings are compatible if they are consistent and equal 
true: 

T{i,a,j,[i) = {Lii V Li2 V ii3)|r.„ A [Lji V L.j2 V ijs)^,,^ = true, 

- where r^y is y-th truth assignment for x-th clause, i.e. y-th strings in the 
truth table for x-th clause. Compatibility box for clauses i and j is the 
following matrix 8x8: 

By = {T{i,a,j,P))sx8- 
The strings' compatibility matrix is the instance: 

G = {Bij)nxm 

- where diagonal boxes are the union matrices of size 8x8. The pattern 
is a box matrix with 8x8 boxes: all elements in the boxes are except 
(1, l)-elements which are 1: 

/ 1 .. 
0.. 



Gi] 



V 



8x8 

Compatibility of system fTTD means that there is a true assignment satisfying 

the given 3-SAT instance. 
2-SAT: The same as the above, except the boxes are 4x4. 
SAT: Let / be a SAT instance: 



f = f\ci, Cj = y 



Let's enumerate literals in each clause. Let's build a compatibility box for 
every two clauses. The box is a rectangular matrix whose elements are 
or 1 depending on whether the appropriate literals in the clauses are 
complimentary: 

Bij = (1 — S{Lia, Ljfs))mxnj 1 

- where 5{a,b) is the Kronecker delta: 

1, a = b 



d{a,b) 



0, a^b 



Compatibility matrix is a box matrix of the compatibility boxes. Each 
clause is presented in this box matrix with one box-row and one box-column 
with the same indexes. The compatibility matrix is the instance: 

G — {Bij)nxn- 

- where diagonal boxes Bu are the union matrices of size rii x Ui. Pattern is 
a box matrix with the same structure but the boxes are filled with except 
their upper left elements which are 1: 

/ 1 

r - 
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Compatibility of system \lM means that the disjunctive normal form of the 

SAT-instance has an implicant. 
The examples are examples of reduction to the Subgraph Isomorphism Problem. 

System 11.11 allows different exact and approximate methods to reduce the num- 
ber of options (the number of non-deterministic solutions to check) or to solve the 
problem. 

For example, let pair (G, S) define a Clique instance. Because matrix multipli- 
cation is a combination of operations "-f" and "x" and all matrices involved in 
system nm are the non-negative matrices, inequalities [LT] will be true for powers of 
matrices G and S, as well. Thus, graph g can be depleted by comparing elements 
of matrix with the value of non-diagonal elements of matrix S^, k > 2. For 
digraphs and k = 2, the procedure can be sketched as follows: 

Step 1: Compute matrix G^ = {aij)„xn', 
Step 2: Rid of all those edges («, j) for which 

Uij < m — 2 V aji < m — 2. 

Iteration of this procedure 0{n — m) times can reduce the instance's dimension or 
even solve the instance. The method can be modified for the Maximum Clique. 
One would start with pattern matrix S of size n x n and then reduce the pattern's 
size using the procedure. 

I. 1. Convex models. Let's fulfill pattern graph s with n ~ m isolated vertices. 
Then, sizes of matrices G and S will be equal and the truncation matrix in formula 

II. 11 will become the union matrix: 



m Pmn = Un- 

Because all matrices involved in system [1.11 are the non- negative matrices and ma- 
trix multiplication involves only summation and multiplication of numbers, the 
system can be rewritten as follows: 

GX > XS 

Let's relax the last system by replacing permutation matrix X with a double sto- 
chastic matrix: 

GX > XS 
< 1 

> 



(1.3) 



Theorem 1.1. NP-instance (G, S) has solution "YES" iff value ^Jn is the solution 
of the following convex program: 



max. 



under constrains\1.3\ 



Proof. Due to the Birkhoff-von Neumann theorem about the double stochastic ma- 
trices [T], value ^/n is the solution of program ll.4l on the Birkhoff polytope 

(1.5) '^^i] = 1) ^^ij 1' ^ij > 0- 
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Matrix X delivers the maximum iff matrix X is a permutation matrix. Permutation 
matrices are the extreme points of the Birkhoff polytope 11.51 □ 



Convex program 11.41 under constrains 11.31 can be solved in polynomial time by 
the ellipsoid method [5] or by the inner point method [12,. 

Convex Drogram I 1 . 41 under constrains ll.Sl is an asymmetric polynomial size model 
for NP-problems. The asymmetry is due to the fixed vertex labeling in which 
adjacency matrix G is written. Relabeling of digraph g will rotate the program's 
solutions over all vertices of the Birkhoff polytope 11.51 

If a NP-problem instance has the negative solution - solution "NO", - then, 
due to theorem II. 1[ there are two cases: the problem's matrix G does not satisfy 
system [T73l at all (an easy case); or, due to the Birkhoff- von Neumann theorem 
about double stochastic matrices and the Carathodory theorem about convex hull, 
the only solutions of the system are certain convex combinations of a permutation 
matrices: 



(1.6) 



2 < a < (n - if + I, 



- because the Birkhoff polytope has dimension {n — if. 

Obviously, constrains 11.31 in theorem 11.11 may be replaced with the following 
system: 

X^G > SX^ 



(1.7) 



Ei < 1 

> 



Theorem 11.11 reduces the Subgraph Isomorphism Problem to the convex pro- 
gramming which is a P-problem [5l I12j . The theorem can be re-formulated: 



Theorem 1.2. NP-instance {G,S) has solution "YES" iff the following convex 
system is compatible: 



GX > XS 

Ej X^j <1 

^ x,j > 

Proof. Value ^/n is the Euclidean norm of the permutation matrices of size nxn. □ 

1.2. Asymmetric linear models. The following example shows that compatibil- 
ity of svstems 11.31 and Fl . 71 is not sufficient for making decisions. 



1 
1 



Exercise 1.3. Let the instance and pattern matrices be as follows: 

^= ( 1 ) ' ^ = 
The only permutation matrices for n = 2 arc 

^1 = ( 1 ^ ) , ^2 = 



1 
1 
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Half-sum of the matrices satisfies systems 11.31 and 11.71 

Xi ^ X2 _ ^Xi ^2 1/1 1 



2 2 '2 2' 2VII 

Nevertheless, instance (G, 5) has solution "NO": value of criterion 1 1.41 on the solu- 
tion is 1 < ^/2. 

The following lemma clarifies the example. 

Lemma 1.4. Let pattern S be a permutation matrix. Let a be the set of all per- 
mutation matrices Y which have the following property: 

Y <G. 

NP-problem instance (G, S) has solution "YES" iff 

3 y e (T : y e ci{S). 



Proof. 



Cl{a) = y Cl{Y). 

Yea 

a 



For example 1, 



G = X2i Cl{S) = {Xi}. 

Let L be any linear functional on i?" . Then, the following asymmetric linear 
program models NP-problems: 

(1.8) L{X) -> max, 

- under constrains 11.31 or 11.71 NP-instance (G, S) has solution "YES" iff there are 
the constrains' extreme points among the optimums of the program. 

If system [TT51 or system [T771 has no solutions at all, then instance (G, S) has solu- 
tion "NO" . If the systems have a solution X and that solution is not a permutation 
matrix, then the solution is a double stochastic matrix. Then, in accordance with 
the Birkhoff - von Neumann theorem: 

X e conv{Xi,X2,...,Xa}, 

- where Xi are permutation matrices and a is in bounds 11.61 The presentation is 
not unique, but it takes O(n^) time to find a permutation matrix which participates 
in one of the presentations: 

Step 1: Select any non-zero element in matrix X. Let Xi-^j-^ be the selection; 
Step 2: Select any non-zero element which is not in row ii nor in column ji. 

Such element exists because X is a double stochastic matrix. Let Xijjj be 

the selection; 

Step 3: Select any non-zero element which is not in rows ii, 12 nor in columns 
ji, j2. Such element exists because X is a double stochastic matrix. Let 
Xi.^j^ be the selection; 

Steps 4~n: And so on until n elements will be selected. 
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Step n + 1: Replace the selected elements with 1 and replace the rest of el- 
ements with 0. The resulting matrix Xi is a permutation matrix which 
participates in the presentation of X as a convex combination of permuta- 
tion matrices. 

If permutation matrix Xi is a solution of system 11.31 or system 11.71 then NP- 
instance (G, S) has solution "YES" . Otherwise, the systems can be fulfilled with 
the following inequality: 

n 

,1=1 

- where and are the indexes selected in the above procedure. In accordance 
with the ellipsoid/separation method [10], the inequality can be used for the next 
iteration of the ellipsoid method . 

1.3. A symmetric linear model. Let's fulfill pattern graph s with n~m isolated 
vertices. Then, sizes of matrices G and S will be equal and the truncation matrix 
in formula ll.il will become the union matrix. Then, system [TTT] can be rewritten: 

(1.9) G>XSX'^ 

- where permutation matrix X is the unknown. 

Let's enumerate all permutation matrices: 

Xl, X2, . . . Xn\- 

Then, integer quadratic system [131 can be replaced with the following integer linear 
system: 

(1.10) 

[ A, e {0, 1} 

- where numbers Xi are the unknown. The system can be relaxed. 
Theorem 1.5. Svstem \l.lO\ is compatible iff the following system is compatible: 

( E.A.X,5Xf <G 
(1.11) { 

[ E.A^ = 1, A, >0 

- where numbers Xi are the unknown. 

Proof. Let numbers solve system [TTTO] and A^^ =1. Then, the following numbers 
are a solution of system [l. Ill 



A, 



1, i = ii 

0, i^ii 

Let numbers A^ solve system [1.111 and A^^ 7^ 0. Because all matrices participated 
in ll.lll are the (0, l)-matrices, 

G ^ Xi , 5* xj" . 



Thus, the following numbers are a solution of system fl. 101 

1, i^ii 
0, iy^ii 



A,; 



□ 
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Due to theorem II .Sj the following linear program will solve NP-problem ll.ll 
(1.12) E^' . ' 

L 

- under constrains ll.lll The program (or its dual, more precisely) can be tried and 
solved with the ellipsoid/separation method [HI [SI [TO]- The separation conditions 
can be arranged with inequalities 

X'lSX'j^ ^ G. 

Due to estimation 1 1.61 the addends on the left side of svstem fl.lll mav be analyzed 
in chunks of size (n — 1)^ + 1. 

Linear program ll.l2^fTTTT] is a symmetric n!-size linear program. The symmetry 
[9] is due to the explicit involvement of all permutation matrices in constrains ll.lll 

For the HC problem, pattern matrix S* is a circular permutation matrix. In this 
case, constrains 11.111 are an explicit expression for the TSP polytope. The ATSP 
with a weight matrix W can be modeled as the following linear program: 

{W, V X^X,SXf) ^ min, 

- under constrains 11.111 The matrix scalar product (*,*) totals products of the 
appropriate elements of its multiplicands. 

According to the Yannakakis theorem [9] , size of the TSP polytope for the sym- 
metric linear program has to be bigger than polynomial. System 11.111 shows that 
n\ is an upper bound for the size. 

1.4. Miscellaneous. 

Lemma 1.6. NP-instance (G, 5*) has solution "NO'' iff the following system is 
incompatible for any permutation matrix R: 

' GX>XS 
X >R 

Xij > 

Proof. If instance (G, S) has solution "YES" , there is permutation matrix X satis- 
fying the system for permutation matrix R = X. Let matrix X satisfy the system 
for some permutation matrix R. X is a double stochastic matrix. Due to the 
Birkhoff-von Neumann theorem, X is a convex combination of several permutation 
matrices. Then, inequality X > R implies that X is a permutation matrix: X = R. 
Then, problem (G, S) has solution "YES" . □ 

Lemma 1.7. Let (G, S*) be such an instance that matrices G and S allow the 
following decomposition: 



G > G1G2, Gi > (0)„x„, G2 > (0) 

^>^1^2, ^1 > (0)„x«, ^2 > (0)„ 



CONVEX AND LINEAR MODELS OF NP-PROBLEMS 



9 



The instance has solution "YES" if the following linear system is compatible: 

' Gi > XSi 

< 1 

x,j > 

Proof. Let matrix X satisfy the system. Then, X is a convex combination of several 
permutation matrices Xf. 

X G conv {Xi \l<i< a}. 

Then, 

G > GiGa > XSiS2X^ > XSX^ > XiSX^ , 

□ 

Any NP-problem (matrix S) has 2" different instances (matrices G). Uhimately, 
if there would exist the mathematical tables of the "YES" /"NO" -instances, the 
tables could be sorted in a way. The binary sorting would reduce the NP-problems 
to the binary search. The computational complexity of the search would be 

0(log2 2"') = 0(n2). 

The "oracle-tables" might be a digital/analog computer which would solve, for 
example, convex program ll.4^1.31 

2. Incidence matrices models 

Let's arbitrarily label/enumerate elements of adjacency matrices G and S: if 
an element is equal a > 0, then that element has a labels (zero-elements have no 
labels). Let's construct in-incidence matrix Iq = {uij): aij = 1 if i-th column of G 
contains j-th label; and cuij = 0, if otherwise. In the same way, let's construct in- 
incidence matrix Ig for matrix S. Also, let's construct out-incidence matrices Oq 
and Os but using rows instead of columns. Direct calculation proves the following 
decompositions: 

(2.1) G = OgII. S = O5/J. 

Let k be the total of all elements of G; and / be the total of all elements of S. 
Then matrices Iq and Oq are nx k; and matrices Is and Os are mx I. In digraph 
terms: numbers k and / are powers of the arc-sets of (multi) digraphs g and s, 
appropriately 0. 

The incidence matrices are total unimodular matrices. They are (0, l)-matrices 
with the following structure: there is one and only one 1 per column. The 0-rows in 
out-incidence matrix indicate sinks in the digraph, and the 0-rows in in-incidence 
matrix indicate sources. The isolated vertices are presented with 0-rows in both 
incidence matrices. From this point of view, the isolated vertices are sinks and 
sources simultaneously. 

In terms of the incidence matrices, svstem fl . II can be rewritten as follows: 

r p^nXOGZp;i = Os 

(2.2) <^ 

I PmnXIaZP^^^Is 



The decomposition can be done for any rectangular matrix. 
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- where permutation matrices X — {xij)nxn and Z — {zij)kxk are the unknown. 
In digraph terms: X presents all n\ ways to label vertices of g, and Z presents all 
k\ ways to label arcs of g. 

Due to the unimodularity of the incidence matrices, each of the two parts of 
system 12.21 has an integral solution if it has a solution at all. The point is the 
existence of such a common integral solution that 

XX^ = Un, ZZ^ - Uk. 

Let's arbitrarily enumerate all permutation matrices of size n x n and k x k: 

Xl, X2, . . . , Xn], Zi, Z2, . . . , Zk\ 

Let's write the following system: 



(2.3) 



Y^i.-i KjPmnXilcZjPi^ — Is 



. A,, e {0, 1} 



where numbers Xij are unknown. 



Lemma 2.1. System \2.2\ has a permutation matrices solution iff system \2.3\ has 



such a solution Xij that 



(2.4) 3 : Ay = 



1, i = ii A j = ji 

0, i^ii V ji 



Proof. Let permutation matrices X and Z he & solution of system [2.21 Then, the 
following numbers are a solution of system [2.31 

1, X,= X A Zj = Z 
0, X,T^ X \J Zj Z 

And yisa-yersa, if a solution of system 12.31 has structure 12.41 then permutation 
matrices Xi^ and Zj^ are a solution of system [521 D 



Ajj 



One might think that solutions 12.41 are the only solutions of system [2731 possible. 
That is incorrect. 

Exercise 2.2. Let instance (G, S) be as follows: 

Og = 



1 






^ 







n 


, Ig = 











,^2 = 1 








I) 


. 


n 



Os - (1), Is = (1) 
The only permutation matrices for n — 2 are 

Xl = 

The following numbers are a solution of system [2T3[ for the instance: 
All — 1, A12 = 1, A21 = 0, A22 — 0. 

Really, 

Pi2XiOgXiP'[2 - (1), P12X1OGX2PI2 - (0) 

P12X1IGX1PI2 = (0), Pi2XiOg^2Pi^2 = (1) 
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Let's add n — m isolated vertices to pattern graph s. That wiU make Pmn = Un 
in systems 12.21 and [2.31 Then, system [2?2] can be rewritten as follows: 

XOgZP£ = Os 



(2.5) 

I V r '7 u 

Ik 

and system [2731 becomes as follows: 



xiGZp;i = Is 



(2.6) 



r E^, \,X,OoZ,Pl = Os 

E.,, x^.xdGZ.Pl = Is 



Ay e {0, 1} 

Lemma 2.3. Any solution of svstem \2.6\ has structure\2.4\ 



Proof. All matrices Oq, Os, Iqi ^-^d Is have one and only one 1 per column. But 
any solution of system [^31 with more than one Ay = 1 will produce, on the left side 
of the system, either a non-(0, l)-matrix or a (0, l)-matrix with more than one 1 
per column. □ 



System 12.61 is a symmetric integer linear model of the NP-problems. The sym- 
metry is due to the explicit involvement of all nlkl combinations of n x n and k x k 
permutation matrices. 

For a given NP-instance, iteration of the following procedure can significantly 
reduce the instance's dimension and even solve the instance in polynomial time: 
Step 1: Build linear combinations of the equations of system [2T6l in order to 

make the right sides of the combinations equal to 0; 
Step 2: Rid system 12.61 of all those Ay which are on the left sides of the 
combinations because, due to lemma [2731 they all are (the ridding can be 
partial when it is difficult to track all Ay = 0). 
From this point of view, the "NO" -instances are such instances for which the result 
of the procedure is that all Ay are 0. 

2.1. A symmetric linear model. System [2761 can be relaxed. 
Theorem 2.4. System \2.6\ has solutions iff the following system has solutions: 

E,,, Ay A,Og^,P;I = Os 



(2.7) 



^vXJgZ.P^, = Is 



EyAy =1, Ay >0 

Proof. Any solution of system 12.61 is a solution of system 12.71 For any solution of 
system 12.71 there is number Xi-^j-^^ > 0. Replacing the number with 1 and the rest 
of numbers Ay with will produce a solution of system l^TBl □ 

The following symmetric linear program models NP-problems: 

^Ay ^ min, 

- under constrains 12.71 In full accordance with the Yannakakis theorem [S] , the 
system has size nlkl. The ellipsoid/separation algorithm [nUTU] can be deployed to 
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solve the program. 

A necessary condition for compatibility of system 12.71 is compatibility of the 
following (nl + fc!)-size system: 

= 1, X, > 
2.2. A convex program. 

Theorem 2.5. Instance (G, S) has solution "YES" iff value ^/n + k is solution of 
the following program: 



(2.8) 

- under constrains 
(2.9) 



max, 



1, E 



1, Zi, > 



Proof. The Birkhoff-von Neumann theorem about the double stochastic matrices 
implies that value y/n + k is solution of program l2. 81 under the following constrains: 



(2.10) 



'E^X^j = 1, E 



E 



1, E 



1, x.i > 



1, z,j > 



The maximum is reachable then and only then when matrices X and Z are permu- 
tation matrices. □ 



Convex program l2.8^2.9l can be solved in polynomial time by the ellipsoid method 
or by the inner point method [3 [12] . 

Linear polynomial size constrains 12.91 are asymmetric. Arc/vertex-relabeling of 
digraph g will rotate vertices of the polytope appropriate to these constrains all 
over the vertices of the polytope defined by equations 12.101 



2.3. Asymmetric linear models. Svstem [2791 alone is insufficient for making de- 
cisions. It misses the following quadratic constrain: 



(2.11) 



ZP^^PikZ^ < Uk. 



Proof. If permutation matrix Z is the "arc" part of a solution of system 12.91 then 
condition l2.11l is true. On the other hand, let matrices Z and X be such a solution 
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of system \2M that inequality 12.111 holds . Let's present double stochastic matrix X 
as a conyex combination of permutation matrices: 

X — Ai^i + X2X2 + . . . , Ai + A2 + . . . = 1, 

- where Xi are permutation matrices. Then, system [2.91 implies : 

G > OcZPi^PikZ^I^ ^ X^OsIgX ^ X^SX > X^SXi. 
Thus, permutation matrix Xj' is a solution of ineauality ll.il □ 

An adequate asymmetric linear model can be build iteratively, as it was done in 
section [TT2] But, let's explore another approach. 

Let's return to symmetric system [^T71 Let cr be a set of all the matrices partici- 
pating on the left side of the system: 

^ = {( xdazM ) M = l>2,. ..,«!; j = 1, 2, . . . , fc!}, 

- where Xi and Zj are permutation matrices of size n x n and k x k, appropriately. 
Let matrix C be in convex hull of the set: 

C e cony(f7) 

For certainty, let matrix C be the center of the polytope: 

r^J_srf x.OgZ.pI \ ^ x,Og z,Pl 

n\k\^\ XdGZjPlk ) nW.y E^XJGE,Z,Pf, 



1,3 



^J_{{{n- l)!)„x„OG((fc - l)!)fexi \ ^ (l)2nxi 

n\k\\iin-iy.)nxnlG{ik-iy.)kxi J n ' 

- because matrices Og and Ig haye one and only one 1 per column. Let r be the 
following set: 

T = {B \ B + C ea}. 

Polytope T is such a shift of polytope a that matrix C gets in the origin of coordi- 
nates. Let matrices Bi be a maximal linear independent subsystem of set r: 

L(i?i,B2,...,i?/5) = L(r) 
^ ' B,eT,l<i<l3 

- where is the linear hull of its arguments. Let us emphasize that 

(2.13) P<2nl. 

Matrices Bi + C, i — 1,2, . . . , (3 are a basis of the minimal hyperplane containing 
set a. 

Theorem 2.6. NP-instance (G, S) has solution "YES" iff the following linear 
system is compatible: 

(2.14) E2/»^'-( )-^' 

- where numbers yi are the unknown. 
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Proof. Necessity. Let instance (G, S) have solution "YES" . Then, system 12.71 is 
compatible. Let's subtract matrix C from both sides of the system; on the left side, 
let's decompose C over numbers A^ ; and let's replace the matrices resulting on the 
left side with their decompositions over basis 12.121 The coefficients resulting on the 
left side are a solution of system [2. 141 

Sufficiency. Let's notice that by the definition of set r, 

V $ e L(t) a^inax >0: fie [0,/imax] ^ ^ conv(r). 

Thus, if system [234] is compatible, then there are such numbers /i and Xij that the 
following system is compatible: 



(2.15) 



xl 



E.,, a.,x,/gZ,f,I = ^,Is + i^Ui 

[ = 1, A., > 



By definition, matrices Os and Is are (0, l)-matrices with one and only one 1 per 
column. Then, the matrix on the right side of system 12.151 has only two different 
elements: 

, ■ 

n n 

By definition, matrices Og and Iq are (0, l)-matrices with one and only one 1 per 
column. Then, all Ay on the left side of system 12.151 will total to two numbers ^ 
and 77, as well: 

" e = (1 - 
7] = /i + (1 - fi)/n 
(n- l)^ + ?7 = 1 
e, ?7 > 0, M > 

Among solutions of the last system, there is the boundary solution (when /i = /imax): 

^ = 1, 77 = 1, e = 0. 

Substitution of = 1 in system [2.151 produces system [^771 Thus, compatibility of 
system 12.141 implies compatibility of system 12.71 i.e. solution "YES" for instance 
iG,S). ^ □ 



Theorem 12 . 141 may be seen as a comparison of given digraph g with the complete 
graph of size n. Due to estimation 12. 131 system [2. 141 is an asymmetric polynomial 
size linear model of NP-problems. The asymmetry is due to the selection of basis 

[m 

The selection of basis 12.1^ is a P-problem. The basis selection is reducible to the 
selection of maximal linear independent subsystems from permutation matrices of 
sizes n X n and k x k, appropriately. The matrices are vertices of the appropriate 
Birkhoff polytopes. Any {n — 1)^ different permutation matrices of size n x n and 
any (fc — 1)^ different permutation matrices of size k x k will produce a polynomial 
size system for the basis selection. The basis selection is reducible to the solution 
of a polynomial size system of linear equations. 
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2.4. Hamiltonian graph. For the Hamiltonian cycle problem [4l[7l[T4j, the pat- 
tern matrix S may be any circular permutation matrix, for example matrix 11.21 
Let's label arcs of cycle s with the indexes of their end- vertices. Then, 

Os = ^, /s = Un- 

Substitution of these matrices in system r2 . 1 41 gives a polynomial size linear system. 
Due to theorem l2.61 the system is compatible iff digraph g is a Hamiltonian digraph. 
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